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Abstract — This  paper  presents  coordination  algorithms  for 
groups  of  mobile  agents  performing  deployment  and  coverage 
tasks.  As  an  important  modeling  constraint,  we  assume  that 
each  mobile  agent  has  a  limited  sensing  or  communication 
radius.  We  focus  on  (1)  a  comprehensive  smoothness  analysis 
of  a  class  of  locational  optimization  functions  (including  a  gen¬ 
eralized  statement  of  the  Conservation-of-Mass  Law),  and  (2)  a 
discrete-time  convergence  result  based  on  a  recently-developed 
generalized  statement  of  LaSalle  Invariance  Principle.  These 
coordination  algorithms  have  convergence  guarantees  and  are 
spatially  distributed  with  respect  to  appropriate  proximity 
graphs.  Numerical  simulations  illustrate  the  results. 

I.  Introduction 

The  current  technological  development  of  relatively  inex¬ 
pensive  communication,  computation,  and  sensing  devices 
has  lead  to  an  intense  research  activity  devoted  to  the 
distributed  control  and  coordination  of  networked  systems. 
The  potential  advantages  of  networked  robotic  systems  are 
their  versatility  and  robustness  in  the  realization  of  multiple 
tasks  such  as  manipulation  in  hazardous  environments,  pol¬ 
lution  detection,  estimation  and  map-building  of  unknown 
environments.  A  fundamental  problem  in  the  distributed 
coordination  of  mobile  robots  is  that  of  providing  stable 
and  decentralized  control  laws  that  are  scalable  with  the 
number  of  network  agents.  Indeed,  there  have  been  various 
efforts  to  provide  rigorous  procedures  with  convergence 
guarantees  using  a  combination  of  potential  energy  shaping 
methods,  gyroscopic  forces,  and  graph  theory  (see  [1],  [2], 
[3],  [4],  [5],  [6],  [7]  and  references  therein).  In  our  previous 
work  [8],  we  studied  distributed  algorithms  for  deployment 
and  optimal  coverage  using  tools  from  computational  ge¬ 
ometry,  nonsmooth  analysis  and  geometric  optimization. 

In  devising  coordination  algorithms  it  is  important  to 
progressively  account  for  the  various  restrictions  that  real- 
world  systems  impose.  Building  on  [8],  this  paper  devel¬ 
ops  spatially-distributed  algorithms  for  coverage  control 
amenable  to  implementation  on  (more)  realistic  network 
models;  we  do  this  by  considering  the  following  new 
aspects.  First,  we  enforce  the  communication  or  sensing 
capacity  of  an  agent  to  be  restricted  to  a  bounded  region, 
typically  much  smaller  than  the  environment  of  interest. 
We  refer  to  these  information  exchanges  between  agents 
as  “limited-range  interactions.”  Second,  we  design  discrete¬ 
time  gradient  ascent  control  laws,  motivated  by  the  argu¬ 
ment  that  discrete-time  feedback  algorithms  are  indeed  the 


ones  truly  amenable  to  implementation  in  a  group  of  agents 
exchanging  information  over  a  communication  network.  We 
deal  with  these  problems  via  a  combination  of  tools  from 
graph  theory,  locational  optimization,  and  systems  theory. 

The  contributions  of  the  paper  are  the  following:  First, 
we  formalize  the  concept  of  spatially-distributed  algorithms 
via  the  notion  of  proximity  graphs.  Second,  we  analyze 
the  smoothness  properties  of  an  important  class  of  objec¬ 
tive  functions,  called  multi-center  functions,  common  in 
locational  optimization,  quantization  theory,  and  geometric 
optimization.  Our  analysis  supersedes  the  results  in  [9], 
[10],  [11],  As  an  important  outcome,  we  determine  the 
extent  in  which  certain  multi-center  functions  are  spatially 
distributed  and  with  respect  to  which  proximity  graphs. 
Finally,  we  design  novel  spatially-distributed  discrete-time 
control  laws  for  groups  of  robots  with  the  objective  of 
steering  the  location  of  a  group  of  robots  to  local  maxima 
of  the  multi-center  function.  We  formally  analyze  their 
performance  and  illustrate  their  behavior  in  simulations. 
In  our  analysis,  we  make  use  of  useful  extensions  of  the 
Conservation-of-Mass  Law  and  of  the  discrete-time  LaSalle 
Invariance  Principle.  These  extensions  are,  to  the  best  of  our 
knowledge,  not  present  in  classical  texts  on  the  subject. 

One  fundamental  scientific  problem  in  the  study  of  coor¬ 
dination  algorithms  is  scalability  with  respect  to  communi¬ 
cation  complexity.  It  is  important  to  design  algorithms  with 
communication  requirements  that  scale  nicely  with  the  num¬ 
ber  of  network  agents.  However,  it  is  impossible  to  quantify 
the  communication  complexity  of  any  algorithm  without  a 
detailed  communication  model.  Adopting  a  computational 
geometric  approach,  this  paper  classifies  the  complexity 
of  coordination  algorithms  via  the  proximity  graphs  with 
respect  to  which  the  algorithms  are  spatially  distributed. 
The  underlying  assumption  is  that  low  complexity  proximity 
graphs  (e.g.,  graphs  with  a  low  number  of  edges)  will  re¬ 
quire  limited  communication  in  a  realistic  implementation. 

Throughout  the  paper  we  consider  only  extremely  simple 
models  for  the  dynamics  of  each  individual  agent  (the  state 
of  each  agent  is  a  point  in  R2  and  its  dynamical  model  is 
an  integrator).  This  feature  is  a  natural  consequence  of  our 
focus  on  network-wide  coordination  aspects.  We  shall  also 
interchangeably  refer  to  agent  as  location  or  point. 

The  organization  of  the  paper  is  as  follows.  In  Section  II 
we  review  various  preliminary  mathematical  concepts  and 
we  develop  useful  extensions  of  classical  analysis  results. 
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In  Section  III  we  study  the  smoothness  of  the  multi-center 
function  and  show  in  what  sense  its  gradient  is  spatially 
distributed.  In  Section  IV  we  design  spatially-distributed 
coverage  algorithms  in  discrete  time.  The  numerical  out¬ 
comes  of  the  algorithms’  implementation  are  reported  in 
Section  V.  Finally,  we  discuss  possible  avenues  of  future 
research  in  Section  VI.  We  refer  the  interested  reader  to  the 
journal  version  of  this  paper  [12]  for  a  detailed  discussion  of 
the  computational  geometric  and  algorithmic  aspects  of  this 
work,  and  the  proofs  of  all  statements  in  this  manuscript. 

II.  Preliminary  developments 

In  this  section  we  present  a  variety  of  preliminary  con¬ 
cepts.  Proximity  graphs  from  computational  geometry  will 
allow  us  to  introduce  the  notion  of  spatially-distributed 
algorithms.  We  refer  the  reader  to  [13,  Chapter  1]  for  an 
exposition  of  standard  graph-theoretical  notions. 


Fig.  1 .  Delaunay,  r-disk,  r-Delaunay,  and  r-limited  Delaunay  graphs  for 
a  random  confi  guration  of  16  generators  with  coordinates  contained  in  the 
rectangle  [0, 1.9]  x  [0,  .75].  The  parameter  r  is  taken  equal  to  .45. 


fixed  V ,  there  exist  n  numbers  M,  (r)  >  0,  *  £  {1, ... ,  n}, 
of  distinct  arcs  arciji(r), . . . ,  arc.j>Mi(r)(r)  of  radius  |  in 
d(Vi(V)  fl  B^{pi))  with  the  property  that 


A.  Voronoi  partitions  and  proximity  graphs 

We  start  by  reviewing  the  notion  of  Voronoi  partition;  we 
refer  the  reader  to  [14],  [10]  for  comprehensive  treatments. 
A  covering  of  R2  is  a  collection  of  subsets  of  R2  whose 
union  is  R2;  a  partition  of  R2  is  a  covering  whose  subsets 
have  disjoint  interiors.  Let  V  be  a  set  of  n  distinct  points 
{pi, . . .  ,pn}  in  R2.  The  Voronoi  partition  of  R2  generated 
by  V  is  the  collection  of  sets  {Vi (’P)}ig{1....,n}  defined  by 
Vi(V)  =  {q  £  R2  |  \\q~pl\\  <  |k -Pi||  ,  for  all  Pj  £  V). 
Here,  ||  •  ||  denotes  the  standard  Euclidean  norm.  It  is 
customary  to  refer  to  Vi(V)  as  Vj. 

Next,  we  present  some  relevant  concepts  on  proximity 
graph  functions  [15].  Let  F(R2)  be  the  collection  of  finite 
point  sets  in  R2;  an  element  of  F(R2)  is  denoted  by 
V  =  {pi , . . .  ,pn}  C  R2,  where  pi, . . .  ,pn  are  distinct.  Let 
G(R2)  be  the  set  of  undirected  graphs  whose  vertex  set  is  an 
element  of  F(R2).  A  proximity  graph  function  Q  :  F(R2)  — > 
G(R2)  associates  to  a  point  set  V  an  undirected  graph  with 
vertex  set  V  and  edge  set  £g{V),  where  £g  :  F(R2)  — > 
F(R2  x  R2)  is  such  that  £giV)  C  V  x  V  \  diag(P  x  V) 
for  any  V.  Here,  diag(P xV)  =  {(p,p)  £VxV  \  p  G  V}. 
In  other  words,  the  edge  set  of  a  proximity  graph  depends 
on  the  location  of  its  vertices.  For  p  £  R2  and  r  £  R+  = 
[0,+oo),  let  Br{p)  =  {q  G  R2  |  \\q  —  p\\  <  r}.  Consider 
the  following  proximity  graphs  (see  Figure  1): 

(i)  the  Delaunay  graph  Qo  with  ( Pi,Pj )  £  £d(P)  if  and 
only  if  Vi(P)  n  Vj(V)  ^  0; 

(ii)  the  r-disk  graph  Gdkk{r)  with  {pi,pj)  £  £disk (P,r)  if 
and  only  if  j|p,;  —  pf  \  <  r; 

(iii)  the  r-Delaunay  graph  ddisknD(f)  is  the  intersection  of 
Sdisk(r)  and  C/D; 

(iv)  the  r-limited  Delaunay  graph  (/ld(t)  with  ( Pi,Pj )  £ 
£ldCP,0  if  and  only  if  A ij(V,r)  =  (VfV)  n 

Remark  2.1:  Let  V  be  a  set  of  n  distinct  points 
{pi, . . .  ,Pn}  C  R2,  and  r  £  R+.  The  collection  {ViiV)  fl 
Br  (pklieli ,...,n}  is  a  partition  of  the  set  U iBr(Pi)  c  R2. 
The  boundary  of  1 f(V)  fl  _EG(p,;)  is  the  union  of  a  finite 
number  of  segments  and  arcs;  see  Figure  2.  Therefore,  at 


d(Vi{V)nBr{Pi))  =  (up.eAAgLD(r)iPi(P)  Ay(r)) 

U  arcj,;(r))  , 

where  A kLD(r),pi  (’^’)  denotes  the  set  of  neighbors  in 
QuriV^r)  of  the  vertex  p,. 


Fig.  2.  The  shaded  regions  are  examples  of  sets  Vi(V)  (T  Br_(pi).  In 
the  fi  gure  on  the  right,  d\{(V)  f  i  B  r  (p):)  does  not  contain  any  arc. 


B.  Spatially-distributed  functions  and  set-valued  maps 

The  notion  of  proximity  graph  is  defined  for  sets  of 
distinct  points  V  =  {pi, . . .  ,pn}.  However,  we  will  often 
consider  tuples  of  elements  of  Bd  of  the  form  P  = 
(pi, . . .  ,pn),  i-e-.  ordered  sets  of  possibly  coincident  points. 
Let  :  (Rd)”  — »  F(Rd)  be  the  natural  immersion,  i.e., 
V  =  iw(P)  is  the  point  set  that  contains  only  the  distinct 
points  in  P.  Note  that  is  invariant  under  permutations  of 
its  arguments  and  that  the  cardinality  of  iF(P)  is  in  general 
less  than  or  equal  to  n.  More  precisely,  if 

S  =  {P  £  (R2)"  |  Pi  =  pj  for  some  if^j  £  {!,... ,  ?r}}, 

then  <  n  if  P  £  S  and  =  n  if  P  qL  S.  The  Voronoi 
covering  V(pr,...,pn)  =  {^(pi,  •  •  ■  ,Pn)}ie{i,...,n}  gen¬ 
erated  by  the  tuple  (pi, . . .  ,pn )  is  defined  by  assigning  to 
each  point  p,  its  corresponding  Voronoi  cell  in  the  Voronoi 
partition  generated  by  V.  Coincident  points  in  the  tuple 
(pi, . . .  ,p„)  have  the  same  Voronoi  cell. 

We  can  now  extend  the  notion  of  proximity  graphs  to 
this  setting.  Given  a  proximity  graph  function  Q  with  edge 
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set  function  £g,  we  define  (with  a  slight  abuse  of  notation) 

g  =  goi¥-.  (r2)"^g(r2), 

£g=£goiv:  (R2)n  -4  F(R2  x  R2). 

Moreover,  the  map  Mg  :  (R2)"  — >  (F(R2))n  is  the 
function  whose  jth  component  is  Mgjipi, . . .  ,pn)  = 
MgtPj(iw(pi,  ■  ■  ■  ,pn))-  Coincident  points  in  the  tuple 
(pi, . . .  ,pn)  have  the  same  set  of  neighbors. 

Given  a  set  Y  and  a  proximity  graph  function  Q,  a  map 
T  :  (R2)n  — »  Yn  is  spatially  distributed  over  g  if  there 
exist  a  map  T  :  R2  x  F(R2)  — >  Y,  with  the  property  that, 
for  all  {pi,...  ,pn)  £  (Rd)”  and  for  all  j  €  {1, . . . ,  n}, 

Tj(pi,  ...,pn)=  f{pj,Mgtj{pi,  .  .  .  ,Pn ))  , 

where  Tj  denotes  the  jth-component  of  T.  In  other 
words,  the  jth  component  of  a  spatially  distributed  map  at 
(pi, . . .  ,pn )  can  be  computed  with  only  the  knowledge  of 
the  vertex  pj  and  the  neighboring  vertices  in  the  undirected 
graph  g{{pi, . . .  ,pn}).  We  are  now  in  a  position  to  state 
an  important  property  of  the  ?’-limited  Delaunay  graph. 

Lemma  2.2:  Let  r  €  R+.  The  set  of  neighbors  map 
A/pLD(r)  of  f?LD  (r)  is  spatially  distributed  over  (/disk(r). 
Loosely  speaking,  this  result  means  that,  to  compute  g ld, 
each  agent  needs  to  know  only  the  location  of  all  the  agents 
in  a  disk  of  radius  r.  This  property  is  to  be  contrasted  with 
the  centralized  computation  required  to  determine  f/disknD- 

C.  Piecewise  smooth  sets  and  Conservation-of-Mass  Law 

A  set  S  C  R2  is  strictly  star-shaped  if  there  exists  a 
point  p  €  S  such  that  for  all  s  €  dS  and  all  A  €  (0, 1], 
one  has  Ap  +  (1  —  A)s  €  int(Sj.  A  curve  C  in  M2  is 
the  image  of  a  map  7  :  [a,b\  — >  R2.  The  map  7  is 
a  parameterization  of  C.  A  curve  7  :  [ a,b\  — >  R2  is 
simple  if  it  is  not  self-intersecting,  i.e.,  if  7  is  injective 
on  (a,  b).  A  curve  is  closed  if  7(a)  =  7 (6).  A  set  O  C 
R2  is  piecewise  smooth  if  its  boundary.  Oil,  is  a  simple 
closed  curve  that  admits  a  continuous  and  piecewise  smooth 
parameterization  71s1  — >  R2.  Likewise,  a  collection  of  sets 
C  R2  |  x  €  (a,  b)J  is  a  piecewise  smooth  family  if 
fl(a;)  is  piecewise  smooth  for  all  x  £  {a,  b),  and  there  exists 
7  :  S1  x  (a,  b)  — >  R2,  (0,x)  1— >  y(0,x),  differentiable  with 
respect  to  its  second  argument  such  that  for  each  x  £  (a,  b), 
the  map  0  1— >  7 x{ff)  =  7 {9,  x)  is  a  continuous  and  piecewise 
smooth  parameterization  of  dQ(x).  We  refer  to  7  as  a 
parameterization  for  the  family  {fl(cc)  C  R2  |  x  £  (a,  6)}. 
The  following  result  is  an  extension  of  the  integral  form 
of  the  Conservation-of-Mass  Law  in  fluid  mechanics  [16], 
Given  a  curve  C  parameterized  by  a  piecewise  smooth  map 
7  :  [a,  b]  ->  C,  recall  that  fcf  =  f* /(7(f))  (7(f)  |  dt 
denotes  the  line  integral  over  C  of  /  :  C  C  R2  — ■>  R, 
and  is  independent  of  the  selected  parameterization. 

Proposition  2.3:  Let  {S7(rc)  C  Q\x  £  (a,  6)}  be  a 
piecewise  smooth  family  with  Yl(x)  strictly  star-shaped  for 
all  x  £  ( a,b ).  Let  <f>  :  Q  x  (a,  b)  — >  R  be  continuously 
differentiable  with  respect  to  its  second  argument  for  all 
x  £  ( a,b )  and  almost  all  q  £  fl(x),  and  such  that  for 
each  x  £  (a,b),  the  maps  q  1— >  tj>{q,x)  and  q  1— >  §4.{q,x) 


are  integrable  on  f l(x).  Then,  the  function  (a,  b)  9  x  1— > 
Jflfxj  x)dq  is  continuously  differentiable  and 

if  </>(<?, X)dq=f  ^{q,x)dq+  f 
ax  JQ(x)  Jn(x)  °x  Jdct(x)  ax 

where  n  :  dfl(x)  — >  R2,  q  1— >  n(q),  is  the  unit  outward 
normal  to  d£l(x)  at  q  £  dfl(x),  and  7  :  S1  x  (a,  b)  — >  R2  is 
a  parameterization  for  the  family  C  R2  |  x  £  {a,  6)}. 

D.  Discrete-time  LaSalle  Invariance  Principle 

We  review  here  some  notions  on  discrete-time  algo¬ 
rithms  [17].  An  algorithm  on  S  C  Rw  is  a  set- valued  map 
T  :  S  — >  2S\{0}.  A  map  from  S'  to  S'  is  simply  a  singleton¬ 
valued  map.  For  any  Xq  £  S,T  generates  feasible  sequences 
as  follows:  given  xn  £  S,  the  map  T  yields  T( xn)  C  S. 
From  this  set,  an  arbitrary  element  xn+i  may  be  selected, 

xn+i  £  T{xn) ,  n  G  N  U  {0}  .  (1) 

An  algorithm  T  is  closed  at  x  £  S  if  for  all  convergent 
sequences  Xk  — ►  x,  x'k  — >  x'  with  x'k  £  T{xk ),  one  has 
x'  £  T{x).  An  algorithm  is  closed  on  W  C  S  if  it  is  closed 
at  x,  for  all  x  £  W.  In  particular,  every  continuous  map 
T  :  S  — >  S  is  closed  on  S.  A  set  C  C  S  is  weakly  positively 
invariant  with  respect  to  T  if  for  any  x$  £  C  there  exists 
x  £  T{x 0)  such  that  x  £  C.  The  function  U  :  S  — >  R 
is  non-increasing  along  T  on  W  if  U{ x')  <  U{x)  for  all 
x'  £  T{x )  and  all  x  £  W. 

Theorem  2.4  (Discrete-time  LaSalle  Invariance  Principle): 
Let  T  be  a  closed  algorithm  on  W  C  Rw  and  let  U  be 
a  continuous  non-increasing  function  along  T  on  W.  Let 
Xq  £  W  and  assume  the  sequence  {xn  \  n  £  N  U  {0}} 
defined  via  a:n+i  £  T(xn )  belongs  to  W  and  is  bounded. 
Then  there  exists  c  £  R  such  that 

xn  — >  M  n  t/-1(c) , 

where  M  is  the  largest  weakly  positively  invariant  set  con¬ 
tained  in  [x  £W  \  3y  £  T{ x)  such  that  U{y)  =  U{x)}. 

III.  Limited-range  locational  optimization 

Let  Q  be  a  simple  convex  polygon  in  R2  including  its 
interior.  Let  diam(Q)  =  maxgiPgQ  ||g  —  p\\.  For  6  <  e  £ 
R+,  let  D[s,e](p)  =  {q  £  R2  |  S  <  ||g  —  p\\  <  e}  denote  the 
annulus  in  R2  centered  at  p  of  radius  <5  and  e;  we  also 
define  D[0<e]{p)  =  Be{p)  =  {q  £  R2  |  ||g-p||  <  e}  and 
D[s,+oo)(p)  =  {<?  e  R2  |  5  <  ||g  -p||}.  Let  nBe<J>){q )  de¬ 
note  the  unit  outward  normal  to  Be(p)  at  q  £  dBe(p).  Given 
S  C  Q.  let  lg  denote  the  indicator  function,  lg(g)  =  1  if 
q  £  S,  and  lg(g)  =  0  if  q  (jL  S.  Throughout  the  paper,  given 
a  point  set  V  =  {pi, . . .  ,pn}.  we  consider  the  restriction  of 
the  Voronoi  partition  V('P)  generated  by  V  to  the  polygon 
Q.  {Vi(V)r\Q}iefi For  ease  of  exposition,  we  denote 
this  restriction  in  the  same  way  as  the  standard  Voronoi 
partition.  Given  a  tuple  P  =  (pi,...,pn)  €E  Qn,  recall 
that  V  =  iw(P)  denotes  the  associated  point  set  that  only 
contains  the  corresponding  distinct  points. 

A  density  function  (f>  :  Q  — »  R+  is  a  bounded  function 
on  Q.  Given  a  set  S  C  Q,  let  area^Sj  =  Js(j){q)dq.  A 


4 


performance  function  f  :  R+  — >  M  is  a  non-increasing  and 
piecewise  differentiable  function  with  finite  jump  disconti¬ 
nuities  at  Ri, ,  R,m  £  R+,  with  Ri  <  •  •  •  <  Rrn .  For 
convenience,  we  set  Ro  =  0  and  Rm+i  =  +oo,  and  write 

m+1 

/(»)  =  £  /“W  l[Ra-i,Ra)ix)  > 

a=l 

where  fa  :  [Ra_i,Ra]  — >  R,  a  G  {1, . . . ,  to  +  1}  are 
non-increasing  and  differentiable  with  fa(Ra)  >  fa+ i(Ra) 
for  a  £  {l,...,m}.  Given  a  density  function  f  and 
a  performance  function  /,  we  consider  the  multi-center 
function  TL  :  Qn  — >  R  defined  by 

n(P)=  [  max  f(\\q-Pi\\)<j>(q)dq.  (2) 

Note  that  Ti  is  an  aggregate  objective  function  since  it 
depends  on  all  the  locations  pi,. . . ,  pn .  It  will  be  of  interest 
to  find  local  maxima  for  TL. 

Remark  3.1:  In  the  optimal  placement  problem  of  large 
numbers  of  spatially-distributed  sensors,  (1)  Ti  provides  the 
expected  value  of  the  sensing  performance  provided  by  the 
group  of  sensors  over  any  point  in  the  environment  Q, 
where  (2)  <j>  is  the  distribution  density  function  representing 
a  measure  of  information  or  probability  that  some  event  take 
place  over  Q,  and  (3)  f  describes  the  sensing  performance 
of  the  sensors.  Because  of  noise  and  loss  of  resolution,  the 
sensing  performance  at  point  q  taken  from  the  /th  sensor  at 
the  position  p*  degrades  with  the  distance  ||g  —  pj||.  • 

Theorem  3.2:  Given  a  density  function  <l>  and  a  perfor¬ 
mance  function  /,  the  multi-center  function  TL  is  globally 
Lipschitz  on  Qn,  and  continuously  differentiable  on  Qn\S. 
where  for  each  i  £  ?x} 


dTL 

dpi 


f)  m 

-^-fih-PiWm^dq  +  ^ifdRa)- 

(P)  dPi  a— 1 


fa+l(Ra))  ’ 


nPRa  ( Pi 
,k(2Ra) 


(3) 


with  arc,:jfe(2f?a),  k  £  {1, . . . ,  M,;(2PQ)}  the  arcs  in  the 
boundary  of  V,  ( P)  0  Bna  (pf). 

This  result  is  a  consequence  of  Proposition  2.3;  we  refer  the 
interested  reader  to  [12]  for  a  detailed  proof.  For  particular 
choices  of  performance  function,  the  corresponding  multi¬ 
center  function  and  its  gradient  have  different  features. 

a)  Centroid  problem:  If  /  is  piecewise  differentiable 
with  no  jump  discontinuities,  then 

W,(P)  =  Iv.mWlmq~Vtm',)dq' 

This  is  the  result  known  in  the  locational  optimization 
literature  [10],  [9],  For  f(x)  =  — x 2,  TL  reads 

n  «  n 

R(p)  =  -Y  /  h-Pi\\2H<i)d<i - Jv*’Pn 

i= l  JvdP)  <=  t 

where  Jw,p  denotes  the  polar  moment  of  inertia  of  the  set 
W  C  Q  about  the  point  p.  In  addition,  the  gradient  of  TL  is 

—  (P)  =  2  f  (q  -  pi)(j)(q)dq  =  2Mv.(p)(CMv.(p)  —pf). 
aPi  JVi(P) 


Here  M w  and  CM  (4-  denote,  respectively,  the  mass  and  the 
center  of  mass  with  respect  to  <j>  of  W  C  Q.  The  critical 
points  of  TL  are  P  £  Qn  with  p,  =  CMy.(P)  for  all  i  £ 
{1, . . .  ,n},  i.e.  centroidal  Voronoi  configurations ,  cf.  [9], 

b)  Area  problem:  If  f(x)  =  l[0,K](a;),  then  TL  corre¬ 
sponds  to  the  area,  measured  according  to  <f>,  covered  by 
the  union  of  the  n  balls  Bp>(pi), . . . ,  Bji(pn),  that  is, 

TL{P)  =  area^ (U™=1  Pr (pi ) ) . 


In  this  case,  the  first  term  in  equation  (3)  vanishes  and  then 


dTL 

dpi 


( P ) 


MR2R) 

/M  /  nBR(Pi) 

fe=l  1  a.ia,ic(2R) 


Given  P  £  Qn,  if  the  /th  agent  is  surrounded  by  neighbors 
in  the  graph  Gld(2R)  in  such  a  way  that  M,  (2f?)  =  0,  then 
the  TL  does  not  depend  on  p, .  This  situation  is  depicted  in 
Figure  2  (see  the  example  on  the  right)  and  captures  the 
fact  that  the  total  area  covered  by  the  agents  is  not  affected 
by  an  infinitesimal  displacement  of  the  zth  agent. 

c)  Mixed  centroid-area  problem:  If  /  is  given  by  a:  1 — >■ 
/ 0)  =  -x2  1[0,r)(x)  +  b  ■  l[Ri+oo)(a;),  for  b  <  - R 2,  the 
multi-center  function  takes  the  form 


n 

W)=-£  >Ap(P)n.BH(pi),pi  +  b  area^Q  \  U "=1Pr(p*))  > 
i— 1 


and  its  partial  derivative  with  respect  p,  is 
dTL 


dp, 


(P)  —  2My.(p)nBR(p.)(CMy.(p)nBR(!)i)  — Pi) 

Mi(2R) 

—  (p2  +  b)  /  nBR(Pi)(t>- 

p=1  •/arci,fc(2fi) 


In  the  particular  case  when  b  =  —R2,  the  function  /  is 
continuous  and  therefore  the  gradient  of  TL  takes  the  form 

dl~L 

—  (P)  =  2My.(p)nBR(p.)(CMy.(p)nBR(j,.)  -Pi)  . 

In  this  case,  the  critical  points  of  TL  are  P  £  Qn  such 
that  pi  =  CM  Vi(P)nBR(pi)  f°r  all  *  G  {1,  We 

refer  to  them  as  R-centroidal  Voronoi  configurations.  For 
R  >  diam(Q),  they  coincide  with  the  standard  centroidal 
Voronoi  configurations  over  Q. 

Corollary  3.3:  Let  <b  and  /  be  a  density  and  a  perfor¬ 
mance  function,  respectively.  The  gradient  of  TL  with  re¬ 
spect  to  the  agents’  location  P  £  Q"  is  spatially  distributed 
over  Gd-  Furthermore,  if  f(x)  =  b  for  all  x  >  R,  then 
the  gradient  of  TL  with  respect  to  the  agents’  location  is 
spatially  distributed  over  Gld( 2.R). 


IV.  Design  of  spatially-distributed  algorithms 

FOR  COVERAGE  CONTROL 

Here,  we  design  discrete-time  implementations  of  the 
gradient  ascent  for  a  general  aggregate  objective  function  TL. 
We  start  by  extending  the  definition  of  TL  to  consider  general 
partitions  W  of  Q.  Let  P  £  Qn  and  let  W  =  {  Wj  C  <Q}”=1 
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be  a  partition  of  Q  such  that  W,  is  a  convex  polygon  and 
Pi  £  Wi,  for  i  £  {1, . . . ,  n}.  Define 

n  ,, 

7fe(P,yV)  =  V  /  f(\\q~Pi\\)<l>(q)dq. 

i= 1  Jw * 

The  function  Pi,  is  differentiable  with  respect  to  its  first 
variable  on  the  whole  Qn.  Note  that,  if  Pi\  (p,W)  = 
fw  /(Ik  ^  p||  )(/>(<?)  d#,  then  we  can  also  write 

n 

He{P,W)  =  '52'H1(pi,Wi). 

i= 1 

The  following  two  equalities  hold 

H(P)  =  He(P ,  V(P)) ,  for  all  P  £  Qn  ,  (4) 
^(P,V(P))=f(P),  for  all  P  £  Qn  \  S  .  (5) 

opi  &Pi 

Let  P  £  S  and  consider  a  partition  W  =  {Wi  C  Q}"=1 
of  Q  such  that  W-L  is  a  convex  polygon  and  p,  £  Wi,  for 
i  £  {1, . . .  ,n}.  Let  i0,  j0  £  {1, . . .  ,n},  i0  ±  j0  such  that 
Pi0  =  Pjo ■  Then,  Vio(P)  =  Vjo(P),  and  V(P)  is  no  longer 
a  partition  of  Q,  but  a  covering.  Nevertheless,  one  could 
consider  the  line  determined  by  the  edge  Wl(1  n  Wj0  and 
the  associated  hyperplanes  HiQtWio nWj(j  and  HjotWiQ nWjo 
such  that  Wio  C  iT,:0,wio nwjo  and  Wj0  C  Hjot wiQnwjo- 
With  a  slight  abuse  of  notation,  redefining 

Ti0(P)  =  VSo(P)  n  Hio,wig nWjo  , 

V3o(p)  =  vjo(p )  n  Hjo,Wi0nwjo  , 

the  collection  V(P)  can  be  seen  a  partition  of  ().  This 
procedure  can  be  extended  if  there  are  more  than  two 
coincident  agents  {*i, . . . ,  is}  at  a  point  p  £  Q  by  defining 

ViM(p)  =  vi„(p)n  (nI/e{i,...,s}\{M}  , 

for  p  £  {1,. In  the  following,  such  a  construction 
will  be  tacitly  performed  for  P  £  S  and  a  partition  W  of  Q. 
One  can  show  fcf.  [12])  that  the  Voronoi  partition  is  optimal 
within  the  set  of  partitions  of  Q,  that  is,  He(P,W)  < 
7 Ye(P,  V(P))  for  all  P  £  Qn.  Moreover,  the  inequality  is 
strict  if  /  is  strictly  decreasing  and  the  partitions  V(P)  and 
VV  differ  by  a  set  of  non-zero  measure. 

We  are  now  ready  to  characterize  a  class  of  algorithms 
with  guaranteed  convergence  to  the  set  of  critical  points  of 
the  aggregate  objective  function  H. 

Proposition  4.1  ( Discrete-time  ascent):  Let  T  :  Qn  — > 
2(i  be  a  closed  algorithm  on  Qn  \  S  such  that: 

(a)  for  all  P  £  Qn,  T{P)  n  5  =  0; 

(b)  for  all  P  £  Qn\S,  P’  £  T{P)  and  i  £  {1, . . . ,  n}, 
W1(p{,Vi(P))>Wi(p<,Vri(P)); 

(c)  for  all  P  G  S  and  P'  £  T(P),  H(P')  >  W(P); 

fd)  if  P  £  Qn  \S  is  not  a  critical  point  of  Pi,  then  for 
all  P'  £  T{P),  there  exists  j  £  {1, . . .  ,n}  such  that 
H1(p'j,Vj(P))>Pi1(pj,Vj(P)). 

Let  Po  £  ()"  denote  the  initial  agents’  location.  Then, 
any  sequence  {Pn  |n€NU  {0}}  generated  according  to 
equation  (1)  converges  to  the  set  of  critical  points  of  Pi. 

This  result  is  a  consequence  of  Theorem  2.4;  we  refer 
the  interested  reader  to  [12]  for  a  detailed  proof.  In  what 


follows,  we  devise  a  general  algorithm  T  :  Qn  — >  2°- 
verifying  properties  (a)-(d)  in  Proposition  4. 1 .  One  can  show 
that  if  po  £  d W,  then  dHlgpn  -  (po)  7^  0  points  toward  the 
interior  of  W.  If  po  £  int(VL)  is  not  a  critical  point,  then 
one  also  has  that  °'HlgpW>  (po)  7^  0.  For  both  cases,  there 
exists  e  =  e(po,  W)  >  0  such  that  the  point  ps  defined  by 


PS  =  Po  +  $ 


dPijjpW) 

dp 


(po)  e  W 


has  the  property  that  Pii(ps)  >  Pi\ (po),  for  all  S  £  (0,  e), 
and  Pii(pt)  =  PL\{po).  As  it  is  usually  done  in  nonlinear 
programming  [17],  the  computation  of  the  step-size  e  can 
be  implemented  numerically  via  a  “line  search”.  With  this 
discussion  in  mind,  let  us  define  the  line  search  algorithm 
T|s  :  Qn  — >  2^"  as  follows: 


Given  P  €  Qn,  let  P'  £  Tls(P)  if,  for  i  £ 
{l,...,n}  with  the  property  that  pi  ^  pj,j  £ 


Pi  =  Pi +  $ 


dH^VijP)) 

dp 


(. Pi )  > 


(6) 


with  S  £ 


e(Pi,Vi(P))  e(Pi,Vi(P)) 
3  ’  2 


,  and,  for  each 


set  {*i, . . . ,  is}  of  coincident  indexes  at  p  £  Q, 


Pi.,  =  Piu  + 


15  ¥P 


(7) 


with  6  £ 


e(Pi.L) 


e(Pi.Yi) 

2 


is  a  partition  of  (P)  =  •  • 

P  G  for  p  £  ,s}- 


where  {Yix , . . . ,  Yis } 
•  =  Vit(P)  verifying 


Proposition  4.2:  The  algorithm  T\s  :  ()"  — >  2(>-  defined 
by  equations  (6)-(7)  is  closed  on  Qn  \  S,  and  verifies 
properties  (a)-(d)  in  Proposition  4.1. 

Corollary  4.3:  Consider  the  algorithm  T\s  :  Qn  — >  2°- 
defined  by  equations  (6)-(7).  Then 


(i)  For  a  general  performance  function  /,  the  algorithm 
Tjs  is  spatially  distributed  over  CD.  If,  in  addition, 
/( x)  =  b  for  all  x  >  R,  then  T[s  is  spatially 
distributed  over  C/ld(2P); 

(ii)  The  sequence  of  agents’  locations  generated  by  2js 
according  to  equation  (1)  starting  at  Pq  £  Qn, 
converges  asymptotically  to  the  set  of  critical  points 
of  the  aggregate  objective  function  Pi. 

Remark  4.4:  In  a  practical  setting,  the  sensing  and/or 
communication  capabilities  of  each  network  agent  are  re¬ 
stricted  to  a  bounded  region  specified  by  a  finite  radius 
r  >  0.  Therefore,  as  a  way  of  approximating  the  gradient 
ascent  of  the  multi-center  function  Pi  for  an  arbitrary  per¬ 
formance  function  /,  one  could  implement  the  discrete-time 
algorithm  7[s  corresponding  to  the  multi-center  function 
induced  by  the  performance  function  fr  :  R+  — >  R  given 
by  fr(x)  =  f(x)  for  x  <  §  and  fn(x)  =  /(diam(Q)) 
for  x  >  | .  This  latter  algorithm  has  the  advantage  of  being 
spatially  distributed  over  C/ldM-  We  do  not  follow  this  path 
here,  and  instead  refer  the  interested  reader  to  [12]  for  a 
complete  exposition,  including  a  detailed  characterization 
of  the  approximation  properties  of  this  procedure.  • 
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V.  Simulations 


VI.  Conclusions  and  future  work 


To  illustrate  the  performance  of  the  coordination  algo¬ 
rithms,  we  include  some  simulation  results.  The  algorithms 
are  implemented  in  Mathematica  as  a  library  of  routines 
and  a  main  program  running  the  simulation.  We  show  exe¬ 
cutions  of  the  discrete-time  algorithm  Tjs  (cf.  equations  (6)- 
(7))  for  the  area  problem  and  the  mixed  centroid-area 
problem  with  discontinuous  sensing  performance  in  Figs.  3 
and  4,  respectively.  Measuring  displacements  in  meters,  we 
consider  the  domain  Q  with  vertices 


We  have  presented  novel  spatially-distributed  discrete¬ 
time  algorithms  for  coordinated  motion  of  groups  of  agents. 
Future  research  lines  include  (1)  distributed  implementation 
of  deterministic  annealing  techniques  (which  promise  to 
overcome  local  maxima),  (2)  visibility-based  algorithms  for 
coverage  in  non-convex  environments,  and  (3)  distributed 
algorithms  for  other  cooperative  behaviors  and  sensing 
tasks,  e.g.,  detection,  estimation,  and  map-building. 
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(2.9325, 1.7),  (2.295, 2.1),  (0.85, 2.3),  (0.17, 1.2)}. 

The  diameter  of  this  domain  is  diam(Q)  =  3.37796.  In  all 
examples,  the  density  function  </>  (represented  by  means  of 
its  contour  plot)  is  the  sum  of  five  Gaussian  functions  of  the 
form  5  exp(6(— (x  —  Center)2  —  (y  —  Wenter)2))  with  centers 
(^center,  2/center)  at  (2,  .25),  (1,2.25),  (1.9,  1.9),  (2.35,1.25) 
and  (.1,  .1).  Measured  with  respect  to  <t>,  the  area  of  the 
domain  is  area^Q)  =  8.61656.  Each  agent  operates  with 
a  finite  sensing/communication  radius  equal  to  r  =  .45. 


Fig.  3.  Area  problem  (with  agent  performance  f(x)  =  1  [q_  ri  (x)): 
discrete-time  algorithm  T\s  for  16  agents  on  a  convex  polygonal  envi¬ 
ronment.  The  left  fi  gure  illustrates  the  gradient  ascent  flow  and  the  right 
ft  gure  represents  the  ft  nal  location  and  Voronoi  partition.  For  each  agent 
i,  the  intersection  V,  D  Br  (p,j )  is  plotted  in  light  gray.  After  18  seconds, 
the  value  of  the  multi-center  function  is  approximately  6.28977. 


m 

Fig.  4.  Mixed  centroid-area  problem  (with  discontinuous  agent  perfor¬ 
mance  f(x)  =  — x 2  l[o,  £)  (tr)  —  diam(Q)2  •  lrr  (a:)):  discrete-time 
algorithm  T\s  for  16  agents  on  a  convex  polygonal  environment.  The  left 
fi  gure  illustrates  the  gradient  ascent  flow  and  the  right  fi  gure  represents 
the  fi  nal  location  and  Voronoi  partition.  For  each  agent  i,  the  intersection 
Vi  n  Bt  (pi)  is  plotted  in  light  gray.  After  13.5  seconds,  the  value  of  the 
multi-center  function  is  approximately  —6.803. 


As  mentioned  in  Remark  4.4,  the  execution  of  the 
coordination  algorithm  in  Figure  4  can  be  regarded  as 
a  limited-range  implementation  of  the  gradient  ascent  of 
the  multi-center  function  TL  corresponding  to  the  agent 
performance  f(x)  =  —a;2;  this  performance  function  does 
not  have  any  range  limitation.  To  compare  both  executions, 
we  computed  the  percentage  error  in  the  value  of  H  at  their 
final  configurations,  which  is  approximately  30.7%.  We  also 
observed  that  the  percentage  error  of  the  performance  of  the 
limited-range  implementation  improves  with  higher  values 
of  the  ratio  r/  diam(Q). 
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